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I give a suggestive theoretical argument against a Gaussian behaviour of 
QED a la Nambu-Jona Lasinio, and show how the numerical data for the 
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1. Introduction 

Quantum Electro-Dynamics (QED) is the most standard and popular 
example of a non-asymptotically free gauge theory. Due to the peculiar 
form of its j3 renormalization group function, as obtained in perturbation 
theory, this model suffers from the Landau pole problem, which means that 
the renormalizcd coupling increases with the energy and diverges at some 
finite energy scale. The only way to construct a consistent theory in the 
infinite cut-off limit is to put the renormalized coupling equal to zero and the 
model is therefore trivial. In spite of that, QED is also the most successful 
quantum field theory from a phenomenological point of view. In fact, at 
energy scales much lower than the cut-off scale, QED describes perfectly the 
electromagnetic interactions. There is no doubt that QED incorporates very 
well the main features of electromagnetic interactions between electrons and 
photons in Nature. 

There are several possible issues to the triviality problem of perturbative 
QED: 

i. Non-perturbative contributions to the P function, which would be 
neghgible at small energies, become relevant at large energy scales and change 
the structure of the /3 function. An ultraviolet stable or non-Gaussian fixed 
point could appear in this way as driving to a non-trivial interacting theory 
in the infinite cut-off limit. This issue includes also the possibility that a 
summation of all the perturbative series be able to change the structure 
of the P function. Therefore "non-perturbative contributions" should be 
understood in a wide sense. 

ii. At large energy scales, and in the context of Unification Models for 
gauge interactions, other physical interactions between particles become rel- 
evant; these new interactions are responsible for the necessary changes in the 
/3 function in order to get a non-trivial model. 

iii. The real world is only described by effective theories, which are inter- 
acting at small energy scales compared with the cut-off scale but are trivial 
if the cut-off is sent to infinity. 

It could be that, after all, the perturbative Landau pole problem re- 
mains, even taking into account all "non-perturbative" contributions to the 
(3 function. However, and independently of this, it is very interesting and 
stimulating, at least from a theoretical point of view, to understand if non- 
asymptotically free gauge theories can be consistently constructed against 
the standard prejudice, and to look for the possible existence of non Gaus- 
sian fixed points in QED. This is however a hard business since, contrary to 
what happens in asymptotically free gauge theories, perturbative solutions 
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of the renormalization group equations cannot be used as a guide in these 
investigations. 

Another source of interest in non-perturbative QED comes from the com- 
posite models for the Higgs sector of the Standard ModeL In fact it has been 
argued that a technigauge model able to strongly couple fermions and to 
produce massive fcrmion-antifcrmion bound states at large energy and with 
large anomalous dimensions (not asymptotically free!) could be on the basis 
of the dynamical mass generation mechanism. 

These arc the main reasons that motivated people 12 years ago [1] to 
start investigations on non-perturbative QED and to search for the existence 
of a non-Gaussian fixed point in this model. An important amount of work 
in this direction has been done in these 12 years. References [1-11] are some 
of the most relevant contributions to the field. 

The existence of a continuous chiral transition at finite inverse gauge cou- 
pling /9, separating a strongly coupled phase where chiral symmetry is broken 
from a weak coupling phase where the symmetry is realized, was found in [1] 
for the non-compact model in the quenched approximation. Another very 
interesting and important contribution to the field was the work by Leung, 
Love and Bardeen [3]: using the continuum formulation and solving the 
Schwingcr-Dyson equations for the fermion self-energy in the quenched lad- 
der approximation, they found large anomalous dimensions near the critical 
point for the composite operator ijjip; this implies a scale dimension of 2 for 
■^■^ and strongly suggests that four-fermion interactions become renormaliz- 
able outside perturbation theory. We should therefore include four-fermion 
interactions in a full analysis of non perturbative QED. 

Reference [4] was the first numerical simulation of the compact model 
with dynamical fermions whereas in [5] it was shown that the compact model 
undergoes a first order chiral phase transition and therefore a continuum limit 
cannot be defined for this model. This point deserves further investigations 
since, as recently found [12], a continuous chiral transition appears in the 
compact model after the inclusion of four fermion interactions. 

Reference [6] contains the results of the first numerical simulation of non- 
compact QED with dynamical fermions. In [7] it was shown that the chiral 
condensate follows a power law behaviour near the critical point, as opposed 
to the essential singularity behaviour or Miransky scaling [2]. 

The DESY-Jiilich group claimed in [8] that they were able to fit the data 
for the chiral condensate with a logarithmically improved scalar mean field 
equation of state. This was a five-parameter fit, one of them being the critical 
couphng. Prom the results of these fits they concluded that QED has a trivial 
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continuum limit. Triviality would be described by a Gaussian fixed point, 
the critical behaviour around it being similar to the one of the model. 
In reference [9], two completely independent and uncorrelated calculations 
were done using different numerical approaches and different operators. The 
extracted critical couplings, which were in perfect agreement, did not support 
however the logarithmically improved scalar mean field fit of the DESY- Jiilich 
group. The data of [9] were compatible, as discussed below, with a Gaussian 
behaviour of QED a la Nambu-Jona Lasinio as well as with a pure power law 
equation of state. 

Actually, the high degree of disagreement between the numerical estima- 
tions of the critical coupling in [8] and [9] is surprising, a disagreement that 
in our more recent calculations in the 14"^ lattice is of about 100 standard 
deviations [13]. Since the numerical data for the chiral condensate in both 
calculations were in good agreement, the only source of disagreement must 
be in the analysis of these data, i.e. in the use of different fitting equations 
in [8] and [9]. The DESY- Jiilich group has recently published new results 
on the logarithmically improved scalar mean field fits of QED with higher 
statistics [14] . In this paper the of the fits is also shown. The best re- 
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ported value is — 7.63 for around 70 degrees of freedom. These numbers 
give an extremely low confidence level for these fits and rule out a Gaussian 
behaviour of QED a la A0^. Indeed, as suggested by Horowitz a few years 
ago [10], if the continuum limit of QED is trivial, triviality would manifest 
itself as in the Nambu-Jona Lasinio model rather than as in the model. 

It is simple to understand this statement if we take it into account that 
QED in the infinite gauge coupling limit and Nambu-Jona Lasinio in the 
infinite four-fermion coupling limit are equivalent. Both models have com- 
posite scalars and pseudo-scalars in the spectrum and, as suggested by Kocic 
and Kogut, and corroborated in a numerical simulation of the Nambu-Jona 
Lasinio model with discrete Z2 symmetry [10], triviality in a theory with 
composite scalars seems not to manifest itself in the same way as in a theory 
with fundamental scalars. This result can be shown for the Nambu-Jona 
Lasinio model in a very simple and elegant way by means of the large N 
expansion. 

Reference [11] contains the results of an analysis of the gauged Nambu- 
Jona Lasinio model, i.e. QED plus four-fermion interactions with continuous 
chiral symmetry. A mean field approximation for the fcrmion field was used 
in [11], whereas fiuctuations of the gauge field were taken into account out- 
side of approximations. The finding of a continuous phase transition line of 
critical points in the gauge-coupling-four-fermion-coupling plane, with crit- 
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ical indices taking non-mean field values, at least near the critical point of 
QED [11], was very suggestive. Notwithstanding that we would like to un- 
derstand the nature of the continuum limit of strongly coupled QED outside 
of approximations. 

I will develop in the next section a suggestive argument against a Gaussian 
behaviour of QED a la Nambu-Jona Lasinio. In section 3 I will discuss some 
recent results of the Frascati-Zaragoza collaboration on the chiral equation 
of state of QED at criticahty, and the last section will be devoted to the 
conchisions. Even if I am not going to discuss it here, I would like to state 
that there are other interesting approaches for the pure gauge compact model 
based on a suppression of monopole contributions [15, 16], the results being 
appealing enough to deserve further investigations. 

2. The Effective Action 

As discussed in the previous section, a Gaussian behaviour of QED as in a 
logarithmically improved scalar mean field theory has been ruled out [9, 13], 
also by the very low confidence level of the fits reported in [14]. This notwith- 
standing, a Gaussian behaviour of QED a la Nambu-Jona Lasinio seems to 
be compatible with the data and it is a difficult task to distinguish this be- 
haviour from a non-Gaussian behaviour. Indeed in both cases the effective S 
exponent which measures the response of the system to a external symme- 
try breaking field at the critical point is less than 3. The effective magnetic 
exponent /Sm, which controls the behaviour of the chiral condensate near the 
critical point, is larger than its mean field value of 0.5 in both cases; the data 
for the susceptibility exponent 7 point clearly to the mean field value 7 = 1 
[17]. To differentiate between logarithmic violations to mean field scaling and 
a pure power law scaling is actually very difficult (see next section). This is 
the reason why, before discussing the numerical results, I would like to give 
an argument based on theoretical grounds, strongly suggesting that QED 
and Nambu-Jona Lasinio are in different universality classes. This argument 
is based on the definition of an effective action for the Nambu-Jona Lasinio 
model and on its comparison with the fermion effective action of QED at 
fixed energy density. 

As we have shown in the past [9, 17], the critical behaviour of QED is 
controlled by a singularity of the fermion effective action, which holds at 
a critical value of the gauge energy density E around 1. I will now show 
how we can define an effective action for the Nambu-Jona Lasinio model, the 
strong four-fermion coupling expansion of which has as first contribution the 
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QED effective action. Before tliat, let me review tlie main features of tlie 
construction of tlie fermion effective action in QED. 

The lattice action of non-compact QED with staggered fermions is 



n,iJ,<v ^ n,iJ, n 

(1) 

where (3 = Ffj_,^ is the naive regularization of the continuum e.m. ten- 

sor, r]ij,{n) are the Kogut-Susskind phases and U^{n) are the compact link 
connections. In the chiral limit this action is invariant under the continuous 
chiral C/(l) transformations 

^„ ^„ e-(-^)"^^ -^- ^„ - ^„e-(-^)-" -^-. (2) 
Action (1) can be written in a compact form as 

= I E KM + ^^i^^ (3) 

n,iJ,<v 

with the following form of the Dirac operator 

A = m7 + iA(C/^(n)) (4) 

A(C/^(n)) being a Hermitian matrix. After integrating out the fermion de- 
grees of freedom the partition function is 

Z = j [rfA^(n)]e"2 ^».''<-^'-'(") det A(m, U^{n)) (5) 

In order to construct the fermion effective action at fixed gauge energy 
density we introduce in the previous expression a trivial identity in the fol- 
lowing way [18]: 



j [dA,{n)]dE5 (±^Y.Fl{n) - e-^^n„<.^M-)detA (6) 

where V in (6) is the lattice volume. If we now introduce the density of states 
N{E) at fixed gauge energy density as 

N{E) = / [dA,{n)]6 E F'^M - e) , (7) 
the partition function (6) can be written as 
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Z = I dEN{E)e-^^f'''{ det A)^, (8) 

with the following definition for the mean value of the fermion determinant 
at fixed E: 

..... _ I[dA,in)]detAS{^j:F^^M-E) 

Contrary to what happens in the compact case, the partition function 
(6), (8) is divergent due to gauge invariance and to the use of non-compact 
variables. These divergences, which are avoidable by fixing the gauge, can 
also be regularized, and they cancel when taking vacuum expectation values. 

Since the pure gauge action is a quadratic form of the gauge fields, the 
density of states N{E) can be analytically computed [9]. At large lattice 
volumes it behaves like 

N{E) = CE^I, (10) 

where C is a divergent constant that contains the divergences due to the 
non-compact integration, and which cancel when computing mean values. 

Now taking into account the fact that the fermion determinant is positive- 
definite, we can introduce the normalized effective action Sff^^{E,m) as 

5gf^(£;,m) = -ilog(detA)^, (11) 

which takes a finite value in the thermodynamical limit. Expressions (8), (10) 
and (11) allow us to write the partition function as a simple one-dimensional 
integral over the gauge energy density: 

Z = J dEe^lio9E-vepE-vs^f^''{E,m) _ (^2) 

The only unknown function in the previous expression is the efi^ective 
action S^^^j:^ {E , m) . It contains all the information on the critical behaviour 
of this model. 

The thermodynamics of this model can be exactly solved in the infinite 
volume limit by the saddle point technique [9]. The mean plaquette energy 
for instance will be given by the solution of the saddle point equation, which 
maximizes the integrand in (12). 

Let me now consider the chiral limit. In this limit, the model has a 
continuous phase transition at some finite value of the gauge coupling, in 
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perfect analogy with magnetic systems. Since first-order phase transitions 
are excluded at small numbers of flavours [19, 20] it is obvious that if the 
effective action Sff^f^{E) is a regular function of E the solutions of the saddle 
point equations will also be regular. No singularities in the plaquette energy 
or specific heat will be found at any finite f3. Therefore in order to get a 
continuous phase transition at finite /? we need a singularity in the effective 
action S^^^^{E) at some finite value of E. This was just what we found by 
numerical computation of S'^^^^{E) [9]. 

The investigation of the main features of this singularity by numerical 
methods is rather difficult. However, what is important for the argument I 
am going to develop here is that this singularity controls the critical behaviour 
of the model and determines its critical exponents. 

The next step in this argument is to define an effective action for the 
Nambu-Jona Lasinio model and to compare it with the QED effective action. 
The lattice action for the Nambu-Jona Lasinio model with continuous chiral 
symmetry and in the chiral limit is 



n,n n,fx 

where G is the four-fermion coupling. This action is also invariant under the 
continuous chiral transformations (2). 

We can bilinearize the previous action with the help of an auxiliary com- 
pact U{1) vector field as follows 

+ $:77^(n)(i/;„e<(")V'n+M - ^n+^e-'''-^^Hr^). (14) 

n,fi 

9^{n) are compact phases taking values in the (— tt, tt) interval. Action (14) 
can be written in a compact form as 

Snjl = i^/^NjLie^)^ (15) 
with the following expression for the Dirac operator: 

AivjL = Ao + (4G')U(C)- (16) 
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Ao is the Dirac operator for a free fermion theory, whereas A(^^) is the QED 
Dirac operator with the gauge field replaced by the auxiliary vector field 

The partition function is the path integral over the auxiliary vector field 
configurations of the determinant of the Dirac operator 



'NJL 



j[dd^]det^r,jL{el). (17) 



Since the fermion determinant in (17) is a periodic function of the compact 
variables 6'^ , we can also use non-compact variables in the definition of Zj^jl- 
The divergences that will appear from the use of non-compact degrees of 
freedom can be regularized; they cancel when computing vacuum expectation 
values, as in non-compact QED. 

The next steps are the same as were used before in the definition of the 
QED effective action. First we define the energy density of the auxiliary 
vector field 9^ just as if it were a gauge field with a "e.m. tensor" : 

KM = e^(n) + e^(n + i,) - elin + u) -e^ (n). (18) 
Inserting now in the partition function (17) a S function 



12^ n,,<u 



and an integral over the energy density E of the auxiliary vector field, we get 
with 

S^ffiE, G) = -1 log ( det A^jl)e, (20) 

and where the definition of the mean value ( )^ over configurations at fixed 
energy density of the vector field is the same as given in (9). 

Again in this case the thermodynamics of the system can be exactly solved 
in the infinite volume limit by the saddle point technique. For instance, the 
mean value of the four-fermion term in the action (13) will be given by 

- ^ dS^ff{E,G) , ^ 

{ '^n'il^n'^n+iJ.'^n+i,, ) = (21) 
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evaluated at the solution E = Eq{G) of the saddle point equation for each 
fixed value of the four-fermion coupling G. 

Doing now a cumulant expansion of ( det /S.njl-i )e inserting in it the 
expression (16) for the Nambu-Jona Lasinio Dirac operator, we get a strong 
four-fermion coupling expansion for the Nambu-Jona Lasinio effective action 

S^/f\E, G) = S%f{E) + S,{E, G). (23) 

The first contribution to this expansion 5'^j^(i5) is just the effective 
action of QED, whereas Si{E, G) in (23) vanishes in the G = oo limit. 

We have now at hand everything we need to finish the argument. As 
discussed before, the QED effective action has a singularity that holds at 
a critical value of the energy E around 1 [9, 17] and which controls the 
critical behaviour of the model. Since the QED effective action is the first 
contribution to the Nambu-Jona Lasinio effective action (23) in the strong 
four fermion coupling expansion of it, the QED singularity will also be a 
singularity of the Nambu-Jona Lasinio effective action. Hence, if by running 
the four-fermion coupling G the solution of the saddle point equations Eq{G) 
for this model reaches the critical energy of QED, the critical behaviour of 
the Nambu-Jona Lasinio model will be controlled also by the same singularity 
that controls the critical behaviour in QED. Both models will be in the same 
universality class even if the behaviour far from the critical point could be 
very different. 

What breaks the previous argument is the assumption that the singularity 
of QED is accessible to the Nambu-Jona Lasinio model. In fact since there 
is no kinetic term for the auxiliary vector field in action (14) and since the 
Nambu-Jona Lasinio effective action is bounded, the solution of the saddle 
point equation is Eq{G) = oo for every value of the four-fermion coupling G. 
The singularity of QED is not accessible to the Nambu-Jona Lasinio model. 
The critical behaviour of the last model must be controlled by a singularity 
in the four-fermion coupling G of its effective action ai E — oo. 

The argument I have developed here has a close resemblance with real 
space renormalization group based arguments. Even if it does not consti- 
tute a rigorous proof, it strongly suggests that both models are in different 
universality classes. 

3. The QED Chiral Equation of State at Criticality 

In this section I will discuss some recent results on the behaviour of the 
chiral order parameter at the critical point of QED [13]. The main ingredi- 
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ents of this calculation are the use of the Microcanonical Fermion Average 
approach (MFA) to simulate dynamical fermions [18], and the computation 
of the chiral transverse susceptibility in the Coulomb phase and in the chiral 
limit, without the help of arbitrary mass extrapolations. 

In the MFA approach the massless Dirac operator is exactly diagonalized 
with the help of a modified Lanczos algorithm. In contrast to what happens 
with other standard algorithms such as Hybrid Monte Carlo, computations 
in the chiral limit are therefore allowed in this approach. We only need to 
take care that the commutation of the chiral and thermodynamical limits 
can be done for the operator we are interested in. In the broken phase, 
characterized by a non-symmetric vacuum under chiral transformations, this 
commutation will give wrong results if applied to operators which are not 
invariant under chiral transformations. However the symmetry of the vacuum 
under chiral transformations allows in principle such a commutation in the 
Coulomb phase. This result can be rigorously demonstrated, for instance for 
the transverse susceptibility. The only ingredient in this demonstration is 
the fact that the susceptibility in the Coulomb phase and in the chiral limit 
is always finite, which implies that the spectral density of eigenvalues of the 
Dirac operator p(A) near the origin behaves like A^, with p > 1. Using this 
result, it is simple to show that the transverse susceptibility is a continuous 
function of the bare fermion mass m at m = 0. 

The computation of the chiral susceptibility in the Coulomb phase and 
in the chiral limit allowed us to get very precise measurements of the critical 
couphngs [13, 17]. This is a crucial point since, as discussed in [9], the 
numerical determination of critical indices is very sensitive to the value of 
the critical coupling. As I mentioned in the introduction, the numerical 
estimations of the critical coupling by the Illinois group and the Frascati- 
Zaragoza group [9] were in very good agreement even if numerical approaches 
and techniques used in these calculations were very different. I would like to 
mention here that we have four independent determinations of the critical 
coupling based on the analysis of the fermion effective action [9] , probability 
distribution function of the chiral order parameter [21], Lee- Yang zeros of the 
partition function and chiral susceptibilities [13, 17]. The results of the four 
calculations were in perfect agreement, but the values extracted from the 
computation of the chiral susceptibility are the best from a statistical point 
of view. The reason for that was the observation that chiral susceptibility 
follows an almost perfect hnear behaviour as a function of the gauge energy 
density in the Coulomb phase, even relatively far from the critical point. The 
very high confidence level of the linear fits for this operator are responsible 
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for the very precise measurements of the critical couplings (see Table I). 

Using the values for the critical couplings reported in Table I, fits of 
the data for the chiral order parameter at the critical point were done with 
the two most reliable equations of state. As previously stated, a Gaussian 
behaviour of strongly coupled QED a la has been ruled out by the 
numerical data [9, 13, 14]. I also discussed how, from theoretical prejudices, 
we expect that, if QED is trivial, triviality would manifest itself as in the 
Nambu-Jona Lasinio model. Even if I have given in the previous section a 
strong argument against a Gaussian behaviour of strongly coupled QED a 
la Nambu-Jona Lasinio, I will show here the results of the fits of the chiral 
order parameter with a Nambu-Jona Lasinio equation of state and with a 
pure power law, the last corresponding to a non-Gaussian behaviour of the 
model. 

The Nambu-Jona Lasinio-like equation of state at the critical point is 
parametrized as follows 

A{ 4,^1: flog + C{ f = m. (24) 

The pure cubic term in this equation fixes the scale of the logarithmic vio- 
lations to scaling. The power of the logarithm is set equal to 1, as follows 
from the expansion of the Nambu-Jona Lasinio model [22]. 

The pure power law equation of state has the following simple form 

B{i)i)Y^m. (25) 

The fits of Figs. 1 and 2 have been done in the mass range 0.01 < 
m < 0.08 and 14^ lattice. In both Figs, we plot the left-hand side of the 
corresponding equation of state against the right-hand side. Very high quality 
fits are obtained in the two follows from the results reported in 

Table II, where I show the ^^-r of the fits in 10^, 12^ and 14^ lattices and for 

' d.o.f. ' 

a number of fiavours running from 1 to 8. The errors in this table take into 
account the error in the determination of the critical couphng (see Table I) . 

The conclusion of this analysis is that the data for the chiral condensate 
at the critical point are not enough sensitive to distinguish between the two 
possibilities. Both equations of state, with the corresponding fitted parame- 
ters, are actually very similar in the mass region explored and even at much 
smaller masses. 

The next step in complexity of the numerical analysis is to test the be- 
haviour of the derivative of the order parameter, i.e. the longitudinal sus- 
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ceptibility. We expect this operator, being a higher order derivative of the 
free energy, to be more sensitive to different equations of state. 

In the case of a Nambu-Jona Lasinio-hke equation of state we get for the 
longitudinal susceptibility xl the following equation 

3m- {iJtlj)x~^ = A{ipijf (26) 
whereas for the pure power law the equation is 

3m- (^V')X"^ = (3-5)5 (^z/')^ (27) 

Equations (26) and (27) have two remarkable things. First, the loga- 
rithmic violations to scaling of equation (24) disappear in the corresponding 
equation for the susceptibility (26). This is very welcome from a numerical 
point of view since it is very difficult to distinguish 5 = 3 with logarithmic 
violations driving to a 5e// < 3 from 6 < 3. Second, if we do a log-log plot of 
log (3m — ( ipip )x~^) against log ( ipip ) we expect a linear behaviour in both 
cases (26), (27). The only difference is the slope of the hnear fit, which should 
be 3 in the Gaussian case and less than 3 for the non-Gaussian one. 

Figure 3 contains the results for the log-log plot in the 12^ lattice. All 

the plotted points in this Fig. arc in the mass range (0.015 0.035). The 

straight line in this Fig. is a linear fit. From the slope of this fit we get that 
the mean field-hke equation of state (Nambu-Jona Lasinio) is also ruled out 
since the extracted value for the critical index is 5 = 1.98(16). 

These results, which arc very appealing, unfortunately need to be im- 
proved in larger lattices before they can be taken as a "proof" of the existence 
of a non- Gaussian fixed point in strongly coupled QED. We have not enough 
statistics to check the previous results in the 14^ lattice. This notwithstand- 
ing, the result 5 — 1.98(16), which is very similar to the quenched value 
reported in [23], is far enough from the Gaussian 5 — 3 value to consider it 
as a strong indication of a non-Gaussian behaviour. 

To finish this section I would like to comment on some features of Fig. 3. 
We have found deviations from the linear behaviour of the function plotted in 
this Fig. in both the left-hand and right-hand corners. The right-hand corner 
corresponds to very small bare fermion masses, and we expect deviations from 
scaling in this region due to finite-size effects. This is the region we should 
explore on larger lattices in order to see if the scaling window enlarges to 
smaller fermion masses. 

The physical origin for the deviations from scaling in the left-hand corner 
is different. In this case we have large bare fermion masses and it is most 
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likely that deviations from scaling come from subleading contributions to the 
dominant scaling in equations (25), (27). 

Another surprising feature, which follows from the analysis of the results 
plotted in Fig. 3, is the value S = 1.98(16). This value is quite different 
from the value S = 2.81 obtained for the same lattice size from the fits of the 
chiral order parameter with the pure power law equation (25). To understand 
this discrepancy we should take into account that the fits of equation (25) 
were done in a much larger mass interval than the one used in Fig. 3. As 
stated before, the chiral order parameter is not very sensitive to the different 
equations of state. Playing with the two free parameters of equation (25), 
we can of course reproduce the data for the chiral condensate in the mass 

interval (0.015 0.035) and 6 ~ 2. The susceptibility fits, as expected, are 

much more sensitive to the form of the equation of state, even if the mass 
interval were these fits work well is significantly reduced. 

4. Summary 

I have done in the Introduction of this paper a short review of some of the 
most relevant contributions to the investigation of non-perturbative QED fi- 
ll]. The relevance of four-fermion interactions in this kind of investigations 
[3] has been emphasized. Unfortunately, and due to technical difficulties, 
results for the full model with four-fermion interactions are available only 
in the mean field approximation for the fermion field [11]. These results 
are however very exciting, since a non-Gaussian behaviour appears in this 
approximation [11]. 

A Gaussian behaviour of QED a la \(f)'^ has been ruled out by the results of 
the Illinois and Frascati-Zaragoza groups [9, 13] as well as by the extremely 
low confidence level of the logarithmically improved scalar mean field fits 
of the DESY-Jiilich group [14]. These results, beside theoretical prejudices 
based on the equivalence between QED and Nambu-Jona Lasinio models in 
the strong coupling limit, tell us that the most reliable realization of triviality 
in QED, assuming it is trivial, is a la Nambu-Jona Lasinio. 

I have given in section 2 a suggestive theoretical argument against a 
Gaussian behaviour of QED a la Nambu-Jona Lasinio. The argument is based 
on the definition of an effective action for the last model and on its comparison 
with the fermion effective action of QED at fixed energy density. I have shown 
how the singularity in the QED effective action, which controls the critical 
behaviour of this model, is not accessible to the Nambu-Jona Lasinio model 
and therefore both models are most likely in different universality classes. 
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The previous expectations are corroborated by the resuhs of the fits of the 
longitudinal susceptibility at the critical point of QED reported in section 3. 
The value of the critical index 5 — 1.98(16) extracted from these fits is quite 
far from its mean field value S = 3. The most reliable scenario is therefore 
that in which strongly coupled QED has a non-Gaussian fixed point and, 
provided that hyperscaling is verified, a non-trivial continuum limit. 
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Figure captions 



Figure 1. 0.97(^V' )^log((:^) + 2.12( )^ against m in the 14^ lattice. 
The sohd hne is a fit with a straight hne of slope equal to 1 and crossing the 
origin. 

Figure 2. 2.39( t/'t/' )^'^^ against m in the 14'* lattice. The solid line is a 
fit with a straight line of slope equal to 1 and crossing the origin. 

Figure 3. log (3m — {i>ip)x~^) against log('0'0) at the critical point of 
the 12^ lattice. The solid line is a linear fit. 
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Table captions 



Table I Critical couplings extracted from the susceptibilities in 10^, 12^ 
and 14'' lattices at Nf = 1, 2, 3, 4, 8. 

Table II Results for the constant B and the critical index 6 (pure power 
law equation (25)) and for the two constants A and C (mean field equation 
(24)) at lattice sizes 10, 12, 14 and Nf = 1, 2, 3, 4, 6, 8. All the fits were done 
in the mass interval (0.01 0.08). 
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